Construction of microscopic model for /-electron systems 
on the basis of j-j coupling scheme 
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We construct a microscopic model for /-electron systems, composed of /-electron hopping, 
Coulomb interaction, and crystalline electric field (CEF) terms. In order to clarify the meaning 
of one /-electron state, here the j-j coupling scheme is considered, since the spin-orbit interaction 
is generally large in /-electron systems. Thus, the /-electron state at each site is labelled by fi, 
namely, the z-component of total angular momentum j. By paying due attention to /-orbital sym- 
metry, the hopping amplitudes between /-electron states are expressed using Slater's integrals. The 
Coulomb interaction terms among the /i-states are written by Slater-Condon or Racah parameters. 
Finally, the CEF terms are obtained from the table of Hutchings. The constructed Hamiltonian 
is regarded as an orbital degenerate Hubbard model, since it includes two pseudo-spin and three 
pseudo-orbital degrees of freedom. For practical purposes, it is further simplified into a couple of 
two-orbital models by discarding one of the three orbitals. One of those simplified models is here 
analyzed using the exact diagonalization method to clarify ground-state properties by evaluating 
several kinds of correlation functions. Especially, the superconducting pair correlation function in 
orbital degenerate systems is carefully calculated based on the concept of off-diagonal long-range 
order. We attempt to discuss a possible relation of the present results with experimental obser- 
vations for recently discovered heavy fermion superconductors CeMIns (M=Ir, Co, and Rh), and 
a comprehensive scenario to understand superconducting and antiferromagnetic tendencies in the 
so-called "115" materials such as CeMIns, UMGas, and PuCoGas from the microscopic viewpoint. 

PACS numbers: 71.27. +a, 71.10.-w, 71.70.Ej, 74.25.-q 



INTRODUCTION 



Since the pioneering discovery of superconductivity in 
CeCu 2 Si 2 p elucidation of the mechanism of unconven- 
tional superconductivity in heavy fermion compounds 
has been one of the central issues both in the experi- 
mental and theoretical research fields of condensed mat- 
ter physics. In particular, further discoveries of supers 
conductivity in uranium compounds such _as UBei3 ,EI 
URu 2 Si 2 j3 UPt 3 ,H UPd 2 Al 3 ,El and UNi 2 Al 3 El have trig- 
gered a rapid increase of vigorous investigations on exotic 
properties of superconductivity in /-electron materials. 

Recently there has occurred a second rush of new dis- 
coveries of superconductivity in /-electron compounds. 
One such is the coexistence of superconductivity and fer- 
romagnetism in UGe 2 HI and URhGe.H After this discov- 
ery, much attention has been attracted to the mechanism 
of unconventional triplet superconductivity in /-electron 
systems, since it is naively believed that triplet pairing 
can coexist in the ferromagnetic phase. A second exam- 
ple is the family of Ce-based heavy -fermion supercon- 
ductors CeMIn 5 (M=Ir, Rh, and Co)IE3y Surprisingly, 
CeCoIns exhibits a superconducting transition temper- 
ature T C =2.3K, which is the highest among those yet 
observed for heavy fermion materials at ambient pres- 
sure. On the other hand, Celrlns shows T C =0.4K, which 
is much less than that of CeCoIn5 . Note that CeRhlns is 
an antiferromagnet with a Neel temperature Tn=3.8K at 
ambient pressure, while under high-pressure it becomes 
superconducting at T C =2.1K. It is interesting that the 



ground state properties of these compounds are easily 
changed by transition metal ions. Third is the supercon- 
ductivity of the filled skutterudite compound PrOs4Sbi 2 
with Tq=jl.85K and a mass enhancement factor as large 
as fiflyja A gain, a relatively high T c is surprising, while 
an interesting point is that the ground state is suggested 
to belong to T 3 , leading to a potential role played by 
quadrupole fluctuations in the mechanism of supercon- 
ductivity. Another amazing result is the discovery of su- 
perconductivity in PuCoGa 5 with T c =1 8. 5K.E3 Although 
it is difficult to synthesize this compound and to measure 
various physical quantities, owing to the inclusion of plu- 
tonium, this is truly high T c and its mechanism is rather 
mysterious. 

Among several remarkable features both of new 
and old /-electron superconductors, one interesting 
and important point is the common observation of 
anisotropic Cooper-pairing, experimentally suggested by 
the power-law behavior of physical quantities in the low- 
temperature region. As is well known, the value of the 
power sensitively depends on the node structure of the 
gap function on the Fermi surface. Thus, by analysing 
carefully the temperature dependence of physical quan- 
tities in experiments, it is possible to deduce the sym- 
metry of the Cooper pairs under the group-theoretical 
restrictions. As a result of this phenomenological analy- 
sis, it has been clarified that the superconducting gap in 
/-electron materials has nodes on the Fermi surface, in- 
dicating that non s-wave Cooper pairing occurs in heavy 
fermion superconductors. 

The appearance of non s-wave pairing itself is easily 



2 



understood as an effect of strong short-range correlation, 
if we recall the fact that the Cooper pair is composed of a 
pair of / electrons with heavy effective mass. Since super- 
conductivity in heavy fermion materials is a second-order 
transition phenomenon, the Ginzburg-Landau theory is 
still applicable to /-electron superconductors if we take 
into account the symmetry of the Cooper pairs as well as 
the crystal structure with the help of group-theoretical 
arguments. Such efforts have been quite useful to pro- 
mote the understanding of the nature of unconventional 
superconductivity without considering the mechanism of 
the Cooper-pair formation.E-3 

As for the microscopic aspects of the mechanism of 
unconventional superconductivity, it has at least been 
confirmed that the strong short-range Coulomb interac- 
tion plays a crucial role. For instance, it brings about 
antiferromagnetic spin fluctuations ,113 which induce sin- 
glet d- wave superconductivity,!^ consistent with the node 
structure suggested by some experiments. It is true that 
in some /-electron materials, anisotropic singlet pairing 
is stabilized, but in UPt3, it has beea experimentally 
confirmed that triplet pairing occurs.liH Quite recently, 
also in UMi#Al 3 , a possibility of triplet pairing has been 
discussed.lia Unfortunately, the triplet pair formation as 
well as the coexistence of superconductivity and ferro- 
magnetism cannot be simply explained as the effect of 
antiferromagnetic spin fluctuations. 

One may naively consider that instead, we can include 
the effect of ferromagnetic spin fluctuations, which was 
proposed for the origin of superfluidity in helium 3. How- 
ever, we then immediately face another problem, i.e., how 
to obtain such fluctuations from the strong short-range 
Coulomb interaction. It may be possible to derive this 
based on a realistic electronic band structure, but we 
do not find that paramagnons are always dominant in 
the spin fluctuation spectrum when we survey the mag- 
netic properties of those materials. In fact, in ruthen- 
ate S^RuCUjij which is experimentally confirmed to be 
a triplet superconductor ,E3 instead we find that the in- 
commensurate antiferromagnetic spin fluctuation is eu-. 
hanced, as observed in neutron scattering experiments.Eil 
Namely, in contrast to our naive expectation, param- 
agnons do not seem to play a central role in the oc- 
currence of triplet superconductivity in the solid state. 
Thus, it is still a puzzling and challenging problem to 
clarify what is the key issue which determines the sym- 
metry of Cooper pairs, singlet or triplet, in /-electron 
superconductors. 

In order to understand this point, it is necessary to 
perform a microscopic analysis based on a model Hamil- 
tonian appropriate for /-electron systems, but unfortu- 
nately, we find such a model to be lacking. It may be 
possible to consider this problem based on the single- 
band Hubbard model, which is believed to describe well 
the properties of strongly correlated electron systems. 
In fact, several basic features of high T c cuprates have 
been understood based on the Hubbard model by ana- 
lyzing it using various kinds of analytical and numeri- 



cal techniques.E3 However, this model is too simplified 
to include differences among materials characterized by 
orbital degrees of freedom, although recently, the active 
role of orbitals has been widely recognized in several ma- 
terials as a means for understanding novel magnetism 
and unconventional superconductivity. For instance, in 
manganites^eg orbitals play a primary role for electronic 
propertiesE3 since orbital ordering has been successfully 
observed in manganites b« means of the resonant X-ray 
scattering measurement .Lj In order to investigate sys- 
tems with active orbital degrees of freedom, it is neces- 
sary to analyze the multi-orbital Hubbard model. Such 
an extension of the model does not bring essential dif- 
ficulty, at least in d-electron materials. For instance, in 
4c? electron systems such as ruthenate Sr2Ru04 intro- 
duced above, it has been found that a three i2g-orbital 
tight-binding model can well reproduce thergxperimental 
results, as can band-structure calculations. c3 

However, if we try to construct such a model also for 
/-electron materials, we immediately find a serious prob- 
lem. Since the spin-orbit interaction is fairly strong in 
/-electron systems, the meaning of a single /-electron 
state is unclear. Even if this point is overcome, we find 
another difficulty in the evaluation of hopping amplitudes 
as well as Coulomb interactions based on a good repre- 
sentation for one /-electron state. In addition, we need 
to include correctly the important effects of crystalline 
electric fields (CEF), which are very often essential for 
some /-electron systems. 

In this paper, we attempt to go beyond all those prob- 
lems through use of the j-j coupling scheme, in order to 
construct an appropriate microscopic Hamiltonian for /- 
electron compounds, including the hopping amplitudes, 
Coulomb interactions, and CEF terms. The Hamiltonian 
obtained is a multi-orbital Hubbard-like model. After 
construction of the model, we exhibit a couple of prac- 
tical models, simplified by discarding one orbital degree 
of freedom. Then, as an typical example, one simplified 
model is analyzed by using exact diagonalization tech- 
niques to evaluate several kinds of correlation functions. 
Especially, the superconducting pair correlation function 
is carefully analyzed based on the concept of off-diagonal 
long-range order. Finally, we will discuss the relation 
between the present theoretical results and some proper- 
ties of CeMIns. We also try to give a possible explana- 
tion for superconducting and antiferromagnetic tenden- 
cies in "115" materials including CeMIns, UMGa5, and 
PuCoGa5, based on our microscopic picture. 

The organization of this paper is as follows. In Sec. 2, 
the construction of the microscopic model Hamiltonian 
based on the j-j coupling scheme is discussed in de- 
tail. The Coulomb interaction terms, CEF terms, and 
the hopping amplitudes are evaluated by paying due at- 
tention to /-orbital symmetry. Two types of simplified 
models are also introduced. In Sec. 3, one of the sim- 
plified models will be analyzed by showing the results of 
spin, charge, orbital, and pair correlation functions. A 
prescription to measure the pair correlation will be also 
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discussed. Finally in Sec. 4, the paper is summarized and 
properties of 115-materials such as CeMIn 5 , UMGa 5 , and 
PuCoGas are discussed. Throughout the paper, we use 
units such that h=k^=l. 

II. MODEL CONSTRUCTION 

In this section we construct a microscopic model 
Hamiltonian for /-electron systems in order to discuss 
unconventional superconductivity and novel magnetism 
from the microscopic viewpoint. For this purpose, the 
model must include correctly the itinerant nature of / 
electrons as well as strong electron correlation and the 
effect of CEF. Namely, the Hamiltonian H should be 
composed of three parts as 

H = -ffkin + lij n t + HqeFi (1) 

where iJkin> -ffint, and -ffcEF are /-electron hopping, 
Coulomb interactions among / electrons, and the CEF 
term, respectively. 

In order to express the Hamiltonian in the form of 
Eq. (|l|), it is requested that the one /-electron state 
should be well-defined by using a label fi. The physical 
meaning of /i will be discussed later in detail, but here 
it is convenient to define a second quantized operator to 
express the state as 

|M>1 = <|0>, (2) 

where |0) denotes the vacuum and a\ is the creation 
operator for / electrons in the /estate at site i. 

By using the second-quantized operator a[ defined 
above, the kinetic term can be written as 

•ffkin = ^i/ a i M a i+azy, (3) 

i,a,/i,f 

where is the overlap integral between the /i- and in- 
states connected by the vector a. The Coulomb interac- 
tion term is given in the form of 

H ^ t = \ I (V:»>v'>V')4 f t a iv a i»' a W> ( 4 ) 

where / is the matrix element for Coulomb interactions. 
Finally, the CEF term is expressed by 

•Hcef = ^ B^ai^aiu, (5) 

where indicates the potential energy for / electrons 
which expresses the effect of the CEF. 

Now our tasks in this section are as follows: (i) Spec- 
ification of the meaning of fi for one /-electron state, 
(ii) Evaluation of the hopping amplitude, Coulomb ma- 
trix elements, and CEF coefficients, (iii) Simplification 
of the model for practical purposes. In the following sub- 
sections, all of these points will be clarified in this order. 



A. j-j coupling scheme 

First let us discuss the meaning of \x in the Hamiltonian 
Eq. (Q). To make it clear, it is instructive to start the 
discussion with the single ion problem using two typical 
approaches, the LS coupling and j-j coupling schemes. 
In general, we consider the /" configuration, where n is 
the number of / electrons included on a localized ion. 
In the LS coupling scheme, first the spin S and angular 
momentum L arc formed by following Hund's rules as 
S=y~ i V_-, Si and L=Y^?_, where Sj and ii are spin and 
angular momentum for i-th /-electron. Note here that 
the Hund's rules are based on the Pauli principle and 
Coulomb interactions among / electrons. After form- 
ing S and L, we include the effect of spin-orbit interac- 
tion \L-S, where A is the spin-orbit coupling. Note that 
A>0 for n<7, while A<0 for n>7. Note also that a good 
quantum number to label such a state is the total an- 
gular momentum J , given by J=L-\-S . Following from 
simple algebra, the ground-state level is characterized by 
J=\L-S\ for n<7, while J=L+S for n>7. 

As is easily understood from the above discussion, the 
LS coupling scheme is valid under the assumption that 
the Hund's rule coupling is much larger than the spin- 
orbit interaction, since S and L arc formed by the Hund's 
rule coupling prior to the inclusion of spin-orbit interac- 
tion. It is considered that this assumption is valid for in- 
sulating compounds with localized / electrons. However, 
when the spin-orbit interaction is not small compared 
with the Hund's rule coupling, especially in actinide ions, 
the above assumption is not always satisfied. In addition, 
if the / electrons become itinerant owing to hybridization 
with the conduction electrons, the effect of Coulomb in- 
teractions would thereby be effectively reduced. In rough 
estimation, the effective size of the Coulomb interaction 
may be as large as the bandwidth of / electrons, lead- 
ing to a violation of the assumption required for the LS 
coupling scheme. 

For /-electron systems in which the spin-orbit interac- 
tion becomes larger than the effective Coulomb interac- 
tions, it is useful to consider the problem in the j-j cou- 
pling scheme. Here we emphasize that the j-j coupling 
scheme is also quite convenient for including many-body 
effects using the standard quantum-field theoretical tech- 
niques, since individual /-electron states can be clearly 
defined, as explained below. First, we include the spin- 
orbit coupling so as to define the state labelled by the 
total angular momentum jj for the i-th electron, given 
by ji—Si+£i. For /-orbitals with £=3, we immediately 
obtain an octet with j=7/2(=3+l/2) and a sextet with 
j =5/2 (=3— 1/2), which are well separated by the spin- 
orbit interaction. Note here that the level for the octet 
is higher than that of the sextet. Then, we take into ac- 
count the effect of Coulomb interactions to accommodate 
n electrons among the sextet and/or octet, leading to the 
ground-state level in the j-j coupling scheme. 

Here, the meaning of /i in our Hamiltonian becomes 
clear. In the case of n<7, [i should be the label to specify 
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the state in the j=5/2 sextet, namely, the z-component of 
the total angular momentum j—5/2 and takes the values 
of /i=-5/2, -3/2, • ■ ■ , 5/2. Note here that for 3<n<7, 
j=7/2 octet is not occupied, since we presume that the 
effect of spin-orbit interaction is larger than that of the 
Hund's rule coupling in the j-j coupling scheme. On 
the other hand, for the case of n>7, /i should be consid- 
ered to specify the state in the j=7/2 octet, since j=5/2 
sextet is fully occupied. Note again that spin-orbit inter- 
action is larger than that of the Hund's rule coupling. In 
this paper, we concentrate only on the case of n<7, es- 
pecially the situation with n—1 and 2, corresponding to 
Ce 3+ and U 4+ ions. Thus, in the following, /i indicates 
the z-component of the total angular momentum which 
specifies the state in the j=5/2 sextet. 

B. Coulomb interaction term 

Now let us consider the ground-state multiplet appear- 
ing in the j-j coupling scheme. For this purpose, first it 
is necessary to define the Coulomb interaction 7 in 77 n t, 
which is the sum of two contributions, written as 



7(/j, v; z/, u') = Knvypi - K^ tfl > 



(6) 



where the former indicates the Coulomb term, while the 
latter denotes the exchange one. It should be noted that 
7 vanishes unless fi+u=fj,'+i'' due to the conservation of 
z-component of total angular momentum. The matrix 



element K, 



is explicitly given by 



KfJ-lfJ-2,f*3fJ.4 ~ y t C MlO"C 'fi 2( j'C 'fj_ 3rT lC \ 



G 17 II G \ /*T\ 

X<,Mi - g'^ 2 ~ 2"ll^ 3 ~~ "2' M4_ 2/' ^' 
where C M(T is the Clebsch-Gordan coefficient, given as 



7/2 - n<7 
7 ' 



(8) 



with (—1) for up (down) real spin. The Coulomb 

matrix element among ^-orbitals is given by 

(mi,m 2 \\m 3 ,m 4 ) 

dridra^Ori)^ ^2)312^3^2)^4 ( r i),(9) 

where 3i2=e 2 /|ri — r 2 \ and <p m (j) is the wavefunction 
for the m state in £—3 orbitals, expressed by the product 
of the radial function and spherical harmonic Y 3m . The 
matrix element (mi, m 2 \ | 7713, 7714) can be written in the 
form of 



(mi,m 2 Urns, mi) 



E 

fc=0 



F k c k (mi 1 m i )c k (m 2 ,mz), (10) 



where the sum on k includes only even valuea-(|fcpO, 2, 4, 
and 6), F k is the Slater-Condon parameterE§EZl includ- 
ing the complex integral of the radial function, and Cfc is 



the Gaunt coefficient ,E2H£j which is tabulated in the stan- 
dard textbooks of quantum mechanics (see, for instance, 
Ref. [§). 

When two electrons are accommodated in the 7 =5/2 
sextet, the allowed values for total angular momentum 
J are 0, 2, and 4 due to the Pauli principle. Thus, the 
Coulomb interaction term should be written in a 15x15 
matrix form. Note that "15" is the sum of the basis num- 
bers for singlet (J=0), quintet (J— 2), and nontet ( J=4). 
As is easily understood, this 15x15 matrix can be de- 
composed into a block-diagonalized form labelled by J 2 , 
including one 3x3 matrix for J z —0, four 2x2 matrices 
for J z =±2 and ±1, and four lxl for J 2 =±4 and ±3. 
We skip the details of tedious calculations for the matrix 
elements and here only summarize the results in the foL 
lowing by using the Racah parameters E k (fc=0,l,2),E3 
which are related to the Slater-Condon parameters F k 
as 



E 



E, = 



E-2 



F 



80 

1225" 



-F l 



12 

441" 



-F 



120 
1225" 



-F l 



1225 

For J z =4 and 3, we obtain 



12 :F*-±F*. 



18 

441' 
1 

44T 



-F 



1(5/2,3/2; 3/2, 5/2) = E - 5E 2 , 



and 



7(5/2, 1/2; 1/2, 5/2) = ^0-5^2, 
respectively. For J z =2 and 1, we obtain 

f 4£ 2 , 



1(3/2, 1/2; 1/2, 3/2) 
7(5/2, -1/2; -1/2, 5/2) 
7(3/2, 1/2; -1/2, 5/2) 



E 
Eq, 



and 



7(3/2,-1/2 
7(5/2,-3/2 
7(3/2,-1/2 



-1/2,3/2) 
-3/2,5/2) 
-3/2,5/2) 



-3V5t; 2 , 



Eo — E 2 , 
E Q + 5E 2 , 
-2VWE 2 , 



Finally, for J z —Q, we obtain 



7(1/2,-1/2 
7(3/2,-3/2 
7(5/2,-5/2 
7(1/2,-1/2 
7(1/2,-1/2 
7(3/2,-3/2 



-1/2,1/2) = E + 2E 2 +E U 
-3/2,3/2) = E -3E 2 +Ei, 
-5/2,5/2) = E Q + 5E 2 +Ei, 
-3/2,3/2) = -Ei-3E 2 , 
-5/2,5/2) = Ei-5E 2 , 
-5/2,5/2) = -Ei. 



Note here the following relations: 

7(/i, z/; z/, //) = I(n',v';v,n), 

and 

7(m, v; v', u') = I(-v, -/i; -//, -v'). 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 
(18) 
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By using these two relations and Eqs. 
obtain all the Coulomb matrix elements 

As a typical example, let us consider the f 2 configu- 
ration. In the j-j coupling schemes, two electrons are 
accommodated in the j=5/2 sextet. When we diagonal- 
ize the 15x15 matrix for Coulomb interaction terms, we 
can easily obtain the eigen energies as Eq—5E2 for the 
J=4 nontet, Eq+9E»Jot the J=2 quintet, and E +3Ei 
for the J=0 singlet E3 Since the Racah parameters are 
all positive, the ground state is specified by J— 4 in the 
j-j coupling scheme. In the LS coupling scheme, on the 
other hand, we obtain the ground-state level as 3 if with 
S=l and L=5 from the Hund's rules. On further inclu- 
sion of the spin-orbit interaction, the ground state be- 
comes characterized by J=4, expressed as 3 H 4 in the 
traditional notation. Note that we are now considering 
a two-electron problem. Thus, if we correctly include 
the effects of Coulomb interactions, it is concluded that 
the same quantum number as that in the LS coupling 
scheme is obtained in the j-j coupling scheme for the 
ground-state multiplet. 

In order to understand the physical meaning of Racah 
parameters, let us consider a simplified Coulomb interac- 
tion term. In the above discussion, the expressions using 
Racah parameters are not convenient since they depend 
on the orbitals in a very complicated manner, although 
they keep the correct symmetry required by group the- 
ory. To clarify their meanings, it may be useful to restart 
the discussion from the following interaction form among 
I = 3 orbitals: 

-ffint = U ^ PimTPiml + U' /2 ^ Pirncr P\m' cr 1 
irn i^^^m^m' 

+ J/% J imafim'cr' hma' fim'a, (19) 

i,a",cr' .m^rn' 

where /i m<T is the annihilation operator for an / electron 
with spin a in the m state of angular momentum I at site 
i and p\ ma = j\ ma j\ma- In this equation, we include only 
three interactions; intra-orbital Coulomb interaction U , 
inter-orbital Coulomb interaction U', and the exchange 
interaction J. Note that the relation U=U'+J holds 
among Coulomb interactions to ensure rotational invari- 
ance in the orbital space. 

By using Clebsch-Gordan coefficients, fi ma with real- 
spin a can be related to as 



we can J— 5/2. Explicitly, J? is written as 



+ ^^' a L-l a iM a L' + l a ip')/ 2 ]: 



(22) 



with cf>± = ^/j(j + l)-p(p±l) =V35/4- M (M±l). 

For two electrons in the j=5/2 sextet, based upon the 
simplified Coulomb interaction term, we can easily ob- 
tain the energy levels as U e s— 5Jh/2 for the J=4 non- 
tet, J7 cff +23J H /2 for the J=2 quintet, and f/ cff +35J H /2 
for the J=0 singlet. When we compare these energy 
levels with the results obtained using Racah parame- 
ters, we understand the correspondence such as Eq^Ucs 
and E2^Jn- Namely, Eo is the effective inter-orbital 
Coulomb interaction, while E2 denotes the Hund's rule 
coupling. Note that E\ does not appear, since it is re- 
lated to the pair-hopping interaction which is not in- 
cluded here. 

Finally, we note the smallncss of Jh, given as Jh= J/ 49. 
The origin of the large reduction factor 1/49 is, in one 
word, due to the neglect of j=7/2 octet. In the Coulomb 
interaction term Eq. (|l9|), the Hund's rule term is sim- 
ply written as —JS 2 . Note the relation S=(g,j — 1)J 
with gj the Lande's g-factor. For j—5/2, we easily ob- 
tain gj=6/7, indicating S= — (1/7)J. Thus, the original 
Hund's rule term is simply rewritten as —(J/49) J 2 . 



C. CEF term 

As is well known, the ground-state multiplet is fur- 
ther split due to the effect of the CEF. In this paper 
we take into account the effect of the CEF by simply 
accommodating several electrons in the one-electron po- 
tential. Namely, £? M „ in Eq. (j|) is expressed by the CEF 
parameters for J=5/2. Since the CEF term depends on 
the crystal structure, it is convenient to write down the 
CEF formula for J =5/ '2 in the typical lattice structure, 
expressed as 



ffcEF = B° 4 (6° 4 + 501), 



for a cubic lattice, 



rrtct nO/nO 

-"CEF — °1 U 2 



(23) 



(24) 



1 3 — am 
< T \I = — a\ 



her/2- 



(20) 



Note here that we consider only the j = 5/2 sextet. The 
Coulomb interaction term for j—5/2 is given by 

H int = U cS J2 n w n W ~ J h^ + (35J H /4)M, (21) 

i/x>/i' 

where 71^=0^0^, A^X^ U eS =U'-J/2, J H =J/49, 
and Ji is the operator for total angular momentum with 



for a tetragonal lattice, and 



rrhcx 
"CEF 



BlOl 



(25) 



for a hexagonal lattice. In each, equation above, O™ is the 
Stevens equivalent operator ,0 composed of J z , J+, and 
J_ . The matrix, formulae have been tabulated in a paper 
by Hutchings.Lj The coefficients B™ are conventionally 
called the CEF parameters, which are, in actuality, deter- 
mined by the fitting of experimental results for physical 
quantities such as magnetic susceptibility. 
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Note here that the above formulae have been obtained 
from the case of J=5/2. In general, the CEF term is 
expressed in matrix form, depending on the value of J; 
for J larger than 5/2, higher terms in 5™ should occur. 
However, as already mentioned above, since in this pa- 
per the effect of the CEF is considered as a one-electron 
potential based on the j-j coupling scheme, it is enough 
to use the CEF term for J=5/2. 

In order to check the validity of our viewpoint regard- 
ing the CEF term, let us see, for instance, the case of 
cubic structure, in which we can easily obtain 



(a) 



£±5/2,±5/2 — 6O-B4, 



±3/2, ±3/2 
±1/2, ±1/2 
±5/2,^3/2 



-18052, 
12052, 



(26) 
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=p3/2,±5/2 



and zero for other /x and v. In the following we consider 
the cases of both f 1 and f 2 configurations. 

When one electron exists in /-orbitals, i.e., in the case 
of the Ce 3+ ion, the ground state is easily obtained by 
the diagonalization of Eq. ( |26| ) . Then we immediately ob- 
tain two eigen energies, E{T 7 )=— 24052 for the T 7 dou- 
blet and E(T 8 ) = 120B2 for the T 8 quartet. Thus, the 
ground state is changed due to the sign of 5°, as shown 
in Fig. 1(a). Qualitatively, the ground state is deduced 
from the positions of ligand anions and the shapes of 
wavef unctions. As in the case of the Ce monopnictide, 
when anions are located along the direction of [1,0,0] from 
the Ce 3+ ion, the Tg level is energetically penalized, since 
the corresponding wavefunctions are expanding along the 
axis directions. On the other hand, as in the case of Ce 
hexaborides with anions located along the [1,1,1] direc- 
tions from the Ce 3+ ions, the T 7 levels are higher due to 
the shape of T 7 wavefunction. 

In order to consider the f 2 configurations such as Pr 4+ 
and U 4+ ions based on the j-j coupling scheme, we put 
two electrons in the potential for the J=5/2 CEF term. 
Here some readers may have a naive question regarding 
its validity. To clarify this point, we analyze the CEF 
level schemes for the f 2 configuration by diagonalizing 
#int+-ffcEF- In Fig. 1(b), we show the eigen energies 
as a function of 5° for a set of typical values of Racah 
parameters, E =5, E±=2, and £2=0.5. For 5°<0, the 
ground state is the T 5 triplet, while the Ti singlet ap- 
pears for 5°>0. At 52=0, the ground-state is the J=4 
nontet. Here it should be noted that the ground-state 
wavefunction consists of a mixture of the J=4, J=2, and 
J=0 states. Note also that J x , not J, is a good quantum 
number to label the state, since the rotational invariance 
in J space is broken by the effect of the CEF for J=5/2. 
However, we confirm that the contribution of the J=4 
component is the largest in the ground state wave func- 
tion. 

In the LS coupling scheme, on the other hand, we need 
to consider a 9x9 matrix form for the J=4 nontet. The 
CEF term for J=4 under cubic symmetry is given by 
H^ F = Bl{6l + hdt)+Bl{dl-2ldt). Note that for the 
case of J=4, there appears a higher order term including 
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FIG. 1: (a) Schematic view of level scheme for n=l un- 
der the cubic CEF. (b) Energy levels as a function of Bl for 
n=2 in the j-j coupling scheme, (c) Energy levels as a func- 
tion of 64 (=605°) for J— 4 nontet under the cubic CEF with 
6 6 (=1260Sg)=0.1. (d) Energy levels as a function of 64 for 
J=4 nontet with 66=— 0.1. 



the CEF coefficients as 



6 . Referring to the paper by Hutchings,! 



B±2,±2 = 
-B±2,=p2 = 

for the J z =±2 sector, 



-660(5° -425°), 
180(552 + 2945°), 



we obtain 



(27) 



5 ±3 ,±3 = -1260(5° + 175°), 

B ThTl = 180(35° + 75°), (28) 

5 ±3 , T i = 5 ±liT 3 = 60^(55° + 635°), 

for the J z =±3 and +1 sectors, and 

5 ±4 .±4 = 840(5° + 65°), 

5 .o = 360(35° -75°), (29) 
B±4,o = 50^4 = 60^(5° - 1265°), 



for the J z =±4 and sectors. 
2x2 and one 3x3 matrices, 
obtain four eigen energies as 



Note that we obtain three 
After diagonalization, wc 



E(Tx) = 1680(52 - 605g) (singlet), 

5(r 3 ) = 240(5^ + 3365°) (doublet), 

5(r 4 ) = 840(52 + 65°) (triplet), 

E(T 5 ) = -120(135° + 2105°) (triplet). 



(30) 



As shown in Fig. 1(c), for 5g>0, the ground state is a 
triplet for 52>0, while it is a Ti singlet for 52<0. On the 
other hand, for 5°<0, the ground state is also a triplet 
for 5°>0 and a Ti singlet for 5°<0, respectively, except 
for the region around 52^0, in which case the doublet 
T3 becomes the ground state, as shown in Fig. 1(d). 

When we compare the LS coupling results with those 
for the j-j coupling, the following two points are noted: 
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(i) Although we commonly observe conversion of the 
ground state between a T 5 triplet and a I\ singlet, the 
correspondence with the sign of B° is reversed between 
the LS and the j-j coupling schemes, (ii) The T3 doublet 
does not appear as the ground state in the j-j coupling 
scheme. The point (i) may not be serious, since it is just 
a problem of parametrization. However, point (ii) looks 
strange at first glance. In order to clarify this point, 
it is useful to change the basis in the j-j coupling. As 
shown in Fig. 1(a), in the one /-electron case, for Bl>0 
(<0), the T 7 doublet is the ground state (excited state), 
while the r§ quartet is the excited state (ground state). 
If we introduce pseudo-spin to distinguish the Kramers 
doublet, T? indicates a single orbital, while denotes 
doubly degenerate orbitals. When we put two electrons 
into the situation shown in Fig. 1(a), owing to the Hund's 
rule coupling we obtain three cases: (i) Two electrons oc- 
cupy T 7 for B%>0. (ii) They occupy T 8 for B^<0. (iii) 
One occupies T 7 and the other does r 8 around B^O. 
After algebraic calculations, we confirm that (i) and (ii) 
correspond to a r x singlet and a T 5 triplet, respectively, 
in agreement with the above result. However, in contrast 
to the naive expectation, (iii) results in a T4 triplet due 
to the Hund's rule coupling. Note that this T4 triplet 
does not become the unique ground state, as shown in 
Fig. 1(b), except at B®=0. In this case, the T3 doublet 
appears only as an excited state, since this T3 doublet 
is found to be composed of two singlet states, each of 
which includes a couple of electrons, one in r 7 and the 
other in Tg. As is easily understood, a local singlet state 
in multi-orbital systems is energetically penalized due to 
the Hund's rule coupling, and thus, the T3 doublet is 
always the excited state, higher than the T4 triplet. 

As long as we consider the j-j coupling scheme for 
the j=5/2 sextet, the T 3 doublet never appears as the 
ground state, which is attributed to ignoring the j=7/2 
octet. However, it is not allowed to include the effect 
of this octet in the present j-j coupling scheme, since 
it just contradicts the original assumption. Namely, if 
the Hund's rule coupling is strong enough so that j=7/2 
cannot be neglected, we must not consider the j-j cou- 
pling schemes, but take the LS coupling scheme from the 
outset. Here we have an interesting conclusion: In the 
j-j coupling scheme, T3 never becomes the ground state. 
Only in the LS coupling scheme with Bg<0 and i?2~0, 
the ground state becomes T3. Except for the disappear- 
ance of the T3 doublet, it is concluded that the CEF 
terms can be treated within the j-j coupling schemes. 



D. Hopping term 

Now we consider the hopping motion of / electrons 
based on the j-j coupling scheme. As discussed in the 
previous subsections, the one /-electron state is labelled 
by fi. Within the second-order perturbation, the hopping 
amplitude i*„ includes two contributions, as shown in 
Fig. 2: (a) The direct hopping of / electrons between the 




i-site i-site 



FIG. 2: Schematic views for (a) direct /-electron hopping 
between neighboring sites along a-direction and (b) hopping 
amplitude i' „ via ligand anions. 

fi-st&te at i site and the ^-statc at i + a site and (b) the 
hopping via ligand anions located along the a-direction 
from the /-ion site. In order to evaluate both contri- 
butions, .basically we employ Slater's two-center integral 
method,t£l as has been done in (i-electron systems. In 
fact, the table for two-center integrals has been extended 
to the case including / orbitals by Takegahara et alxLA 
However, as emphasized in the previous subsections, it is 
necessary to pay due attention to the fact that spin and 
orbital are mixed. Thus, in this case, it is not allowed to 
simply use the Slater-Koster table. 

To avoid any confusion, it is convenient to step back 
to /-electron operators in the £=3 multiplet, defined as 
fima- Using Clebsch-Gordan coefficients, f\ ma with real- 
spin tr is related to as 

a ifi = J]] CW./i,Ai-cr/2,£r- (31) 
cr 

Note again here that a is defined to be +1 (—1) for up 
(down) spin. 

First, let us consider the direct hopping of / electrons 
along the a-direction, as shown in Fig. 2(a). Since the 
real spin should be conserved in the hopping process, 
is given as 

= ^WCW^ 1 3 a ^-CT/2;3,t'-CT/2' (32) 

where Tf m . t i , is the hopping amplitude of electrons 
between (£,m)- and (£', m')-states along the a-direction. 

Now the problem is reduced to the evaluation of 
Lf m . t i m i. Although we can simply consult the paper 
of Slater and Koster,0 a convenient formula has been 
obtained by Sharma for the prerlap integral between two 
orbitals, (£, m) and (£',m')E3 connected by unit vector 
a. It is expressed as 

I ~Ajr / 4tt 
Z?m;*',m' = W*) \j \j ^+1^ ^ <P)Yt> m >W, <p), 

(33) 

where (££'a) denotes Slater's two-center integral through 
the a bond, for instance, it is (//tr) for £=£'=3 and (/per) 
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for l=2> and £'—1. Yp m is the spherical harmonic function 
with m=—£, ■■■,£. 9 and ip are polar and azimuth angles, 
respectively, to specify the vector a as 



a = (sin 9 cos <p, sin 9 sin ip, cos ( 



(34) 



Let us consider the direct hopping between / orbitals 
in nearest neighbor sites by putting £=£'=3. After some 
calculations, we obtain the hopping amplitudes as fol- 
lows. For diagonal elements, we obtain 



±5/2,±5/2 
a 

±3/2, ±3/2 



= 5to sin 9, 

= t sin 2 9(1 + 15 cos 2 9), 



, , ., - 2t (l-2cos 2 fl + 5cos 4 6>), 
where the energy unit to is given by 



(35) 



E. Simplified models 

In previous subsections, we have completed the con- 
struction of a model Hamiltonian which is expected to 
be a basic model to investigate the microscopic proper- 
ties of /-electron systems. However, it includes six states 
per site, i.e., three Kramers doublets. For practical pur- 
poses, it is more convenient to simplify the model fur- 
ther by discarding one of the three Kramers doublets. 
Although the hopping amplitudes are obtained in a gen- 
eral form applicable to any type of crystal structure, it 
is necessary to set the crystal structure for the purpose 
of actual calculations. 



1. T 8 model 



to 



!(/M 



(36) 



where (ffcr) is the Slater-Koster two-center integral be- 
tween adjacent / orbitals. Note that t* t _ = 0. For 
off-diagonal elements, we obtain 



±5/2, ±1/2 
a 

±5/2,=p3/2 
a 

±l/2,=p3/2 



= -tox/lOe^ 2 ^ sin 2 6>(1 -3 cos 2 (9), 



hs/Se^v sin 4 9, 
-t \/2e* 2i v sin 2 9(1 



(37) 



5 cos 2 9), 



and 



= t \/l0e- 3 ^ sin 2 6 sin 20, 



''5/2,-1/2 ~~ t l/2,-5/2 

*B/2,3/2 = -i a 3/2 ,_5/2 = - 2< oV5e-^sin 2 0sin26\ 
t* /2 3/2 = -t a 3/2 _ 1/2 = t V2e^ sin 29(1 - 5 cos 2 t 



(38) 

Note that iS M =tJ£. 

Next we consider the hopping of / electrons t'* v 
through ligand anions, as shown in Fig. 2(b). Assume 
that i?-orbitals are situated along the a-direction and /- 
electrons hop to the neighboring /-site through these £- 
orbitals. Again we note that the real spin is conserved 
in the hopping process. In second-order perturbation, we 
can easily obtain 



J -3, f i-<r/2;e,m J -e,m;3M-a/2 
Si -El 



(39) 



where eg and Sf are the energy levels for I- and /-orbitals 
and m runs between —I and I. Due to the orthogonality 
of spherical harmonics, it is easily reduced to 



t' 



t 



(IV 



£f - £i (ffcr) ' 



(40) 



It is stressed that the a-dependence of is the same 
as that of t* although the energy unit is changed from 
(ffcr) to (f£a) 2 / (ef — eg)- Note that in general, the sign 
is also changed. 



First let us consider a simple cubic lattice composed of 
ions with / electrons. Under the cubic CEF, as discussed 
above, for i?°<0 the quartet becomes the ground 
state, while the doublet is the excited state. This 
situation is considered to be realized when ligand an- 
ions are situated in the face-centered or body-centered 
position, due to a naive discussion of the point-charge 
picture. Note that the effect of ligand anions is included 
only through the CEF term, and we do not consider its 
effect on /-electron hopping. Then, it is possible to con- 
sider a model consisting of Ts, neglecting under the 
condition of i?4<0. 

To consider the model including T 8 levels, it is useful 
to define new operators with "orbital" degrees of freedom 
to distinguish two Kramers doublets included in r 8 as 



(41) 



Ci«T = VW6ai-5/2 + V^76 a i3/2 ) 
dai = V^/6a i5/2 + V^76 a i-3/2, 



for a orbital electrons and 



■i/3T 



Oi-1/2) c i/31 — Oil/2, 



(42) 



for {3 orbital electrons, respectively. Note that the T 7 
state is given as a 7 orbital as 



c i7T = V^76«i-5/2 - V^76 a i3/2, 

c ni = V^76 a i5/2 - V576ai_ 3 /2. 



(43) 



For the standard time reversal operator JC=— ia y K, 
where K denotes an operator to take the complex conju- 
gate, we can easily show the relation 



/Cci- 



CTCi. 



(44) 



Note that this has the same definition for real spin. 

For hopping in the xy plane and along the z-axis, we 
simply set (9, ip) to be (tt/2,0), (n/2,n/2), and (0,0) 
for a=x=[l,0,0], y=[0,l,0], and z=[0,0,l], respectively. 
Then, by using the general results in the previous subsec- 
tion, we easily obtain t* v between neighboring / orbitals 
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in the xy plane and along the z axis. Further we trans- 
form the basis by the above definitions for operators with 
orbital degrees of freedom. The results are given asc*3e!J 



t , = t 



for the x-direction, 



1 -1/V3 
-1/V3 1/3 




(45) 



1 1/V3 
t y TT , = t \ 1/V3 1/3 




(46) 



for the y direction, and 




e TT , = t | 4/3 I , 




(47) 



for the z direction. Note that t=6t =(9/28)(ffa). Here 
we note two points: (i) Hopping amplitudes among r 8 
orbitals are iust the same as those for the e g orbitals of 
3d electrons.cil This is quite natural if we recall the fact 
that r 8 is isomorphic to T 3 x Tg, where T 3 indicates E 
representation for the orbital part and Tq denotes the 
spin part, (ii) The r 7 orbital is localized, since the cor- 
responding wavefunction has nodes along the axis direc- 
tion. Since Tr is the excited state for i?4<0 and it is also 
localized, electrons do not occupy the r 7 orbital. 

After lengthy and tedious calculations to transform the 
basis of the Couloeab interaction terms, the Hamiltonian 
for r 8 is given byu3 

H * = trT'cLafi+ar'* -sj^ipia- Pib)/2 

i,a,(7,r,T / i 

+ U PirWir j + U' ^ Piapib 



it a' C\ T > o 



(48) 



where p\ T <j= c\ T<J Ci T(J and PiT=^2 a P\<jt- Note that in the 
second term, the level splitting e is simply introduced 
to include partially the effect of a tetragonal CEF. In 
the Coulomb interaction terms, U, U', J, and J' denote 
intra-orbital, inter-orbital, exchange, and pair-hopping 
interactions, respectively, expressed by the Racah param- 
eters E k as 



U = E + E 1 +2E 2 , 

U' = E a + {2/3)E 2 , 

J = 5E 2 , 

J' = E 1 -(ll/3)E 2 . 



(49) 



Note that the relation U=U'+J+J' holds, ensuring ro- 
tational invariance in pseudo-orbital space for the inter- 
action part. 



We believe that this Tg Hamiltonian provides a simple, 
but non-trivial model to consider superconductivity and 
magnetism in /-electron systems. Note again that it is 
essentially the samepas the model for e g electron systems 
such as manganitesjiil although the coupling with Jahn- 
Teller distortion is not included in the present model. 
Due to the complex interplay and competition among 
charge, spin, and orbital degrees of freedomj-a rich phase 
diagram has been obtained for manganites.c^l Thus, it is 
definitely expected that a similar richness will also be 
unveiled for /-electron systems based on the r 8 model 
Hamiltonian. 



2. T 7 model 

The r 8 Hamiltonian iJ 8 is believed to be a non- 
trivial model for /-electron systems, since we can im- 
mediately grasp some essential aspects of /-electron ma- 
terials viewed as charge-spin-orbital complex systems. In 
fact, based on i? 8 , we can have insight into the effect on 
superconductivity of cooperation and competition among 
spin and orbital fluctuations in the new family of heavy 
fermion materials CeMIns (M=Ir, Cp-and Rh), with the 
HoCoGas-type tetragonal structured 

Among CeMIns materials, the de Haas- van Alphen 
(dHvA) effect-h|as been successfully observed in Celrlns 
and CeCoInsJlj'aboth of which have huge electronic spe- 
cific heat coefficients of several hundreds of mJ/K 2 -mol. 
The angular dependence of major experimental dHvA 
frequency branches is well, explained by a quasi two- 
dimensional Fermi surface .ejuj This is a clear advan- 
tage when we construct a model Hamiltonian, since it 
is enough to consider the two-dimensional case. 

In order to analyze further the electronic properties 
from a quantitative viewpoint, it is desirable to include 
more precisely the effect of the actual crystal structure. 
For instance, in the previous subsection, we have intro- 
duced the level splitting e between T 8 levels in a simple 
way, in order to take into account partially the effect of a 
tetragonal CEF. However, in an actual tetragonal situa- 
tion, there occur three Kramers doublets, two T7 and one 
r 6 . Thus, in this subsection, we attempt to incorporate 
the effect of a tetragonal CEF more faithfully into the 
model. 

First let us consider the tetragonal CEF term. After 
the diagonalization, we obtain three Kramers doublets, 
two and one r 6 , expressed as 



|lf;) = ( P al ±5/2 +qal T3/2 )\0), 
(~ial ±5/2 +pal T3/2 )\0), 



f>) 

|r 6 > 



(50) 



a\ ±1/2 \0) 



The coefficients p and q are determined by the CEF pa- 
rameters as 



V(l + a)/2, q = - «)/2, 



(51) 
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FIG. 3: Schematic view of a two-dimensional lattice com- 
posed of Ce and In ions for CeMIns (M=Ir, Rh, and Co). 



where Q=B\j\B\\ and 



= 12(£° + 20B°)/e 7 , 



with 



e 7 = 12 J(B° + 2QBlf + 20(B1) 2 . 



The corresponding eigen energies are given by 



£(rf) 

^(4 2) ) 

E(T 6 ) 



4(B% - 15^2 
4(5° - 
-8(B§ 



15S°) 



155°) 



£7/2, 



(52) 



(53) 



(54) 



respectively. 

If the r 6 level is higher than the other two r 7 levels, it 
is possible to obtain another simplified model for the r 7 
levels. In fact, in some experimental papers on CeMIns, 
we have found that the position of the-To is the highest 
among the three Kramers doublets. E3o Thus, a simpli- 
fied and realistic model for CeMIns can be constructed 
by including only two T 7 levels. 

In this crystal structure, the Ce ions form a two- 
dimensional tetragonal lattice, while the In ions are lo- 
cated in the face-centered positions of this Ce-ion net- 
work, as shown in Fig. 3. We include the nearest- 
neighbor hopping between Ce 4/ orbitals as well as next- 
nearest neighbor hopping via In bp orbitals. Due to the 
symmetry, the latter process contributes only along the 
diagonal direction, while hopping via p-orbitals along the 
bond direction vanishes. 

Based on the same two-dimensional lattice, we have 
further constructed a Hamiltonian called thp-j'/-p model" 
by keeping the f-p hybridization explicitly^ In order to 
compare the tight-binding dispersion directly with the 
results of a relativistic band-structure calculation, it is 
important to keep the hybridization term. However, here 
we consider an /-electron model wherein we include the 
effects of f-p hybridization in the next-nearest neighbor 
hopping of / electrons. We believe it will be adequate for 



the treatment of superconductivity and magnetism if we 
further add the Coulomb interaction terms to the model. 

Since two r 7 's are given by the linear combinations 
of /i=±3/2 and ±5/2, we obtain the T 7 model just by 
discarding the /i=±l/2 states. It is convenient to define 
additional new operators as 



QlT = a 5/2, Qll = a-5/2, 
Q2T = 0-3/2, c -\2{ = a 3 / 2 - 

Again, we can easily show the relation JCci Ta =aci 
For the nearest neighbor hopping, we obtain 

5 y/E' 
s/E 1 



trr> = (t/6) 



(55) 



(56) 



Note that t has the same definition as for the T 8 model. 
On the other hand, for the next-nearest neighbor hopping 
via the p orbitals, we obtain 



t' TT , = (t'/6) 



5 -V5 

-VE 1 



(57) 



for both [1,1,0] and [1,-1,0] directions. The next- 
nearest neighbor hopping amplitude is given by 



t' 



9 {fpof 
28 Sf - Er 



(58) 



where (/per) is the Slater-Koster two-center integral for 
/ and p orbitals. Note that in general t'<0. 
The r 7 model is then written as 

H 7 = t TT 'c\ Ta Ci +a . T > a + Y *rr' c iTa C i+br' CT 

i,a,cr,r.r' i,b, cr,r,r' 



+ (£7/2) Y ( aTz + - a 2 T x ) TT >c\ Tt7 C iT > a 

\,(T,T.T' 

+ ^ TJrPir^Pirl +U')] PiaPib 
i,r i 

+ J^2pilaPi2a + J' ^ C irT C ir'| CiT i CiT 'T' ( 59 ) 
i.er i,T^r' 

where a and b, respectively, denote the vectors con- 
necting nearest- and next-nearest neighbor sites. t x 
and t z are Pauli matrices. The energy level shift has 
been neglected here for simplicity. The intra-orbital 
Coulomb interactions are given by Ui=Eq— ZE2+E1 and 
U2=Eq+5E2+Ei for the 1 and 2 orbitals, respectively, 
and the inter-orbital Coulomb interaction is denoted by 
U'=Eo—5E2- The Hund's coupling and the pair-hopping 
interactions are given by J=10E 2 and J'=Ei, respec- 
tively. 

Three comments on H-j are in order. (1) There is 
no transverse component in the Hund's coupling term. 
(2) Uiy^U2 since the 1 and 2 orbitals belong to differ- 
ent symmetries. (3) For the case of n=l, both J and 
J' are sometimes neglected for simplicity, since they be- 
come irrelevant in the large- V limit. When we are in- 
terested in the strongly correlated region, it is possible 
to safely set S 1 =£ , 2 =0 for the case of n=l, leading to 
Ui=U 2 =U'=E . 
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III. ANALYSIS AND RESULTS 

In the previous section, we have discussed in detail 
the construction of a microscopic model for /-electron 
systems based on the j-j coupling scheme. For practi- 
cal purposes, we have also prepared a couple of simpli- 
fied models corresponding to an actual crystal structure. 
Those models can be analyzed with several techniques 
from both an analytical and a numerical viewpoint. It is 
impossible to show here the results for all models. In this 
section, we show exact diagonalization results for Hs as 
a typical example of an orbital degenerate system. 

A. Measured quantities 

Before proceeding to the exhibition of calculated re- 
sults, let us discuss the quantities we would measure in 
a model including orbital degree of freedom. First of all, 
in order to see the spin and charge structure, it is quite 
natural to measure the spin and charge correlation func- 
tions, defined as 

5(q) = ^ e ^( i -J)( SziSzj ), (60) 
ij 

and 

C(q)=J2e i «< i -»(p iP5 ), (61) 

ij 

respectively, where s z i — ^2 T (piT^ — /°iri)/2 and p\ = 
53tct Pira- On the other hand, in order to see the or- 
bital structure, it is quite useful to calculate the orbital 
correlation function, given by 

oMq) = I> 4q - (i - j Wjr')- (62) 

i j 

As is easily understood from the above equations, even 
if the orbital degree of freedom is further included in the 
model Hamiltonian, there is no ambiguity in the defini- 
tions for spin and charge correlations, given by sums in 
terms of an orbital index. For the orbital correlations, we 
also consider as an orbital-dependent charge correlation, 
which can be defined without any difficulty. 

However, when we attempt to discuss the supercon- 
ducting pair correlation, we immediately encounter some 
serious problems. Namely, it is necessary to determine 
the form factors for the pairing states in both momentum 
and orbital space. Note that even in the single-orbital 
Hubbard model, it is necessary to determine the momen- 
tum dependence, but in general, our task can be reduced 
through the help of a group-theoretical analysis. In the 
orbital degenerate model, the extra task of determining 
the form factor in orbital space is also imposed upon us. 
Especially in the present realistic models, in which the 
orbital degree of freedom does not provide a good quan- 
tum number, the situation seems to be very severe at 



first glance, since we cannot resort to a group-theoretical 
analysis for the orbital dependence of the form factor. 

In order to overcome this problem, we carefully recon- 
sider it,based on the concept of off-diagonal long-range 
order £3 Let us define the pair correlation as 

P ± = (*±*L>. (63) 

where the subscript — (+) denotes spin singlet (triplet) 
pairing and $± is a pair operator, given in general by 

$± = Yl ^M'( C ^T C i+aM'i ±C iMi C i+aM / T )/\/2. (64) 
i,a,/i,/V 

Since a pair with zero total momentum is considered here, 
p , depends only on a, which is the vector connecting 

the sites. Note also that the orbital dependence of 
is not solely determined by symmetry considerations. 

We can easily show that the pair correlation is written 
in the form 

p± = ^'^'.to'^'. (65) 

a,b,/.i,i/,/i' ,u' 

where the matrix element p^,,,.> k,„,> is written as 

+ ( c iMl c i+aM'T c j+bi/T C jf 

± (ciMlCi + a^Tc]+b,'ii T )]/2. (66) 

Here it is necessaw to recall the concept of off-diagonal 
long-range order .c3 Namely, the occurrence of supercon- 
ductivity is detected when the largest eigenvalue of the 
pair correlation matrix becomes the order of N, where N 
is the number of electrons. Thus, in a small cluster cal- 
culation, the possible superconducting pair state should 
be defined by the eigenstate with the maximum eigen- 
value p m ax of the matrix P. Of course, in order to prove 
the existence of off-diagonal long-range order, we need to 
show that Pmax becomes the order of N with increasing 
the cluster size. However, at some cluster size, a possi- 
ble pairing state is determined without ambiguity by the 
eigenstate with p max . 

A prescription to define the pair operator is briefly 
summarized. First, we determine the pair correlation 
matrix pi,„/ \, vv i- Then we diagonalize it, and the form 
factor <^a„„/ will be defined by the eigenstate correspond- 
ing to the maximum eigenvalue of the pair correlation 
matrix. 



B. Symmetry of the orbital dependent Cooper-pair 

In the previous subsection, we have discussed a way 
to define the pair operator in orbital degenerate systems. 
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Here we further consider the symmetry of the Cooper pair 
in orbital space from a general point of view. Anomalous 
Green's functions for spin triplet and singlet pairs are, 
respectively, defined by 

F^(k,r) = -[(cktfWc^i) + (c k/4 (r)c_ kl/T )], (67) 
and 

^( k ' T ) = -[(cVtCOc-k^) ~ (ck^(-r)c_ kl/T )], (68) 

where c^ a (T)—e Hl o &ili<y er HT and r denotes imaginary 
time. Due to cyclic invariance of the trace, we obtain the 
following relations: 



FU-k,r) = F* (k,/3-r), 



(69) 



where (3 is the inverse temperature. Further, after Fourier 
transformation, we obtain 



F%{~Kiu n ) = F^{k,-iuj n ), 
FU-k,iuj n ) = -FUk,-iu n ] 



(70) 



where uj n is the fermion Matsubara frequency, given by 
uj n =TiT{2n+l) with n an integer. These are the basic 
relations representing the symmetry of Cooper pairs in 
orbital degenerate systems. 

Except for differences in the spin states such as sin- 
glet and triplet, "parity" is the only extra good quantum 
number to specify the pairing state, since there is no con- 
servation law regarding orbital degrees of freedom. Thus, 
we can produce even- and odd-parity states, satisfying 



Pf£(k,iw n ) = ±f£(k,iw n ), 



(71) 



where V indicate the operator to change k to — k. It 
is easy to understand that F± is given by linear combi- 
nations of F^(k,iuj n ) and F^ v (k,—iuo n ). Although we 
cannot define F± uniquely, it is convenient to consider 
combinations of odd and even functions both for orbital 
exchange and sign-change in frequency. If we introduce 
two orbitals "a" and "b" , then after some simple algebra, 
we obtain 



FUKiUn) = J2 A t F t( k ' i "n), 
z=l 



(72) 



where Af is an arbitrary number and F^ 1 (k, iui n ) is given 
by 

F±(k,iu, n ) = F a s a (k,z^„)±F a s a (k,-^„), 
F±{k,iw n ) = F b s b (k,z Wn )±F b s b (k,^„), 



F±(k,iu n ) = F a s b (k,* W „)±F a s b (k,- J c„) 

+ F^(k,ioj n )±F° a (k,-iu n ), 

F±(k,iw n ) = ^ b (k, !U „)TF a s b (k,-H) 

- F b s a (k, JW „)±F b s a (k,^„). 



(73) 



Since the present Hamiltonian is not invariant with or- 
bital exchange^ the terms "orbital symmetric" and 



"orbital antisymmetric" cannot be used as a label to 
specify the state. However, it is convenient to catego- 
rize F: (k, iu> n ) as follows: The i^ + (k, iu) n ) (i=l~3) are 
orbital-symmetric and even function of frequency, while 
the F~(k,iuj n ) (i=l~3) are orbital-symmetric and odd 
function of frequency. On the other hand, F~£ (k, iui n ) is 
orbital-antisymmetric and an odd function of frequency, 
while i^(k, iu> n ) is orbital-antisymmetric and an even 
function of frequency. 

If we consider the zero- frequency limit, we can sim- 
ply drop the odd functions of frequency. Namely, the 
even-parity pair is expressed as a combination of orbital- 
symmetric functions, while the odd-parity solution is 
just given by F 4 ~(k, iu> n ), the orbital-antisymmetric part. 
Note, however, that this is just an accident due to the 
vanishing of odd-functions of frequency. We must not 
consider that the orbital state provides a good label to 
specify the pairing symmetry, although it is a convenient 
way to classify the pair state. Again, it is stressed that 
spin state and parity are labels used to specify the pairing 
symmetry. 

For the case of spin triplet, it is easy to repeat the 
same discussion as that for spin singlet. Note that in 
the zero- frequency limit, the odd-parity solution is given 
by the combination of orbital-symmetric parts, while the 
even-parity function is expressed only by the orbital- 
antisymmetric part. 

Up to now we have discussed the general situation in 
which the system is not invariant under the exchange of 
orbitals. We believe that this is quite natural since in 
actual materials, the shapes of orbitals are different from 
one another. However, it is possible to consider a fic- 
titious Hamiltonian which is invariant under orbital ex- 
change. In such a situation, the words "orbital symmet- 
ric" and "orbital antisymmetric" recover physical mean- 
ings which distinguish the pairing states. Namely, we can 
further classify anomalous Green's functions as 



+8 (Mw„) 
+a (k, iu) n ) 



F 

K. 
Fl s (k,icj n ) 



E-=iA + F+(k,^„), 

F 4 + (k,io;„), 
F 4 (k,iw„), 



(74) 



where the subscripts + and — denotes parity even and 
odd, while "s" and "a" indicate the orbital-symmetric 
and orbital antisymmetric, respectively. In this case, the 
frequency dependence is automatically classified: Even- 
function for +S and —a and odd-function for -|-a and — s. 

Finally, it is again stressed that when we consider pair 
correlation in the static limit, at first glance, orbital sym- 
metric and antisymmetric parts seem to be separated, 
but such a separation is spurious owing to the speciality 
of the static limit. 



C. Numerical Results 

Next we show results obtained using exact diagonal- 
ization techniques for an 8-site tilted cluster. Since we 
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(b) charge correlation 
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FIG. 4: (a) Spin, (b) charge, and (c) orbital correlation func- 
tions vs. e for U=U'=4 and J=J'=0. (d) Ground-state phase 
diagram for the T$ model for U=U' with J=J'=0. (e) Pair 
correlations vs. e both for spin singlet and triplet states, (f) 
Pair amplitude in momentum space for the singlet state hav- 
ing the maximum eigenvalue of the pair correlation matrix. 
Note that it is shown separately for aa, bb, and ab+ba pair 
states. 



attempt to make some comparison with CeMIns , we con- 
sider a two-dimensional lattice. First we consider the 
case of quarter filling, i.e., n=l. Note here that J and 
J' are simply neglected, since they are considered to be 
irrelevant in the case of n=l. As shown in Fig. 4(a), for 
e<0.5, there is no dominant component in S(q), indicat- 
ing a paramagnetic (PM) phase, while for e>0.5, S(ir,ir) 
becomes abruptly dominant, strongly suggesting an an- 
tiferromagnetic (AF) phase. This abrupt change is due 
to level crossing, but it is instructive to see the charge 
and orbital correlations. In Fig. 4(b), we show charge 
correlation, exhibiting again an abrupt change around 
e«0.5. A decrease of charge correlation is the signature 
of a metal-insulator transition, since the charge fluctua- 
tions are suppressed in the Mott insulator. As observed 
in Fig. 4(c), we can see an orbital disordered (OD) phase 
for e<0.5 and a ferro-orbital (FO) phase for e>0.5. By 
monitoring the change in spin correlations, we can draw 
the ground-state phase diagram for n=l, as shown in 
Fig. 4(d). We see OD/PM phase for small U and FO/AF 
phase for large U region. 

In Fig. 4(e), we show the eigenvalues of the pair cor- 



relation matrices for spin singlet and triplet states as a 
function of e. The state with the largest eignevalue is 
a singlet, while a triplet appears as the second largest 
eigenvalue state. Even if we introduce a finite value of J, 
unfortunctely, this triplet state is not the largest eigen- 
value state, at least within the present 8-site calculation. 
However, in the previous analysis using the random phase 
approximation, the appearance ota triplet state has been 
suggested in the small e regionlH3 It may be related to 
the present result having the triplet state as the second 
largest eigenvalue state. 

In Fig. 4(f), the pair amplitudes for the largest eigen- 
value state is shown in momentum space. Since we are 
measuring the static pair correlation, we can observe 
here an "orbital symmetric" state with "aa", "bb", and 
"ab+ba", while the "orbital antisymmetric" part such 
as "ab— ba" is not mixed. As discussed in the previous 
subsection, it is just an accident due to the static limit, 
and the orbital symmetric state is not separated in gen- 
eral. Next we focus on the momentum dependence of 
each component. As seen in Fig. 4(f), the "aa" and "bb" 
parts of the pair amplitude exhibit -Big, while "ab+ba" 
seems to belong to Ai g . If we recall that "a" and "b" each 
correspond to e g orbitals, then they have different signs 
for 7r/4 rotation in space. Thus, there exists a discrep- 
ancy in the symmetry between ( "aa" , "bb" ) and "ab+ba" 
components, although the total symmetry should be Bi g 
in this case. 

We have obtained a simultaneous onset of spin and or- 
bital ordering, which can be interpreted as a transition 
between PM and AFM phases controlled by the suppres- 
sion of orbital fluctuations. Moreover, just around the 
transition regime between PM and AFM phases, we can 
expect the appearance of d-wave superconductivity in- 
duced by AFM spin fluctuations, which has been identi- 
fied as the Bi g pairing state with the largest eigenvalue. 
In fact, the present model has been analyzed using the 
random phase approximation, leading to a transition be- 
tween PM arid d-wave superconducting phases with an 
increase of eo If we further increase e, eventually the 
system becomes an AFM insulating phase, and this suc- 
cessive transition agrees with the present results. Here it 
is stressed that the orbital degree of freedom plays a key 
role to control the ground-state properties. We believe 
that this viewpoint is applicable to the appearance of d- 
wave superconductivity in CeMIns, although the quanti- 
tative features should be discussed in detail using the T7 
model. 

Next we consider the case of n=1.5 for the same F$ 
model. For e=0, we depict the phase diagram in the 
(U',J) plane, as shown in Fig. 5(a). Note that the re- 
gion with J>U' is unphysical and here we just ignore it. 
We see the appearance of the ferromagnetic (FM) phase, 
which is stabilized due to the strong Hund's rule coupling. 
The interesting point is the pair amplitude inside the FM 
phase, as shown in Fig. 5(b). The momentum depen- 
dence of the largest eigenvalue state is characterized by 
even parity. As naively expected, due to the Hund's rule 



14 




FIG. 5: (a) Phase diagram for e=0 and n=1.5. Note that 

the region with U' < J is unphysical. (b) Amplitude for pIG g. Configuration for /-electrons in Ce, U, and Pu ions 
triplet pairing in momentum space, evaluated from the eigen- fof the pseudo gpin repr esentations. Each horizontal bar in- 
state corresponding to the maximum eigenvalue of the pair dicateg the pgeudo Qrbital 
correlation matrix. 



coupling, a local triplet is formed, but in order to gain 
kinetic energy, this pair should be spatially extended, as 
shown in Fig. 5(b). Note here that this is Ai g -like, since 
it is an even-parity state. Since we are now considering 
the pairing state in the strong coupling region, there is 
no reason to exclude the triplet state with even parity, 
which is not directly connected to the usual odd-parity 
triplet state. This is not surprising, since local triplet 
pair must have even parity, provided that one ion with 
/-orbitals is included in the unit cell. Thus, quite gener- 
ally, we can conclude that the odd-parity triplet pair is 
not stabilized in the local Hund's rule coupling, at least in 
the BravaisJattice, which is consistent with the previous 
discussion. E3 

As for the mechanism of the attractive interaction in 
the pair inside the FM phase, it should be considered to 
be orbital fluctuations, since spin fluctuations are dead 
in the fully spin-polarized state. This point is intriguing 
due to a possible connection with the coexistence of su- 
perconductivity and ferromagnetism in UGe2- Of course, 
in order to confirm such an exotic type of superconduc- 
tivity, it would be necessary to consider the problem in 
the bulk limit using, for instance, the Green's function 
method. This case will be discussed elsewhere. 

Finally, we note that due care should be paid to the 
correspondence between pseudo-spin and total angular 
momentum J, when we intend to make a comparison 
with actual materials. In experiments on /-electron com- 
pounds, for instance, the external magnetic field couples 
to J, not to the pseudo-spin. Thus, when we attempt 
to calculate physical quantities observed in actuality, it 
is necessary to consider the response of J. The present 
analysis based on the pseudo-spin model is useful to see 
qualitiave tendencies, while a quantiative discussion is 
left for future development. 



IV. SUMMARY AND DISCUSSION 

We have constructed a microscopic model Hamiltonian 
for /-electron systems based on the j-j coupling scheme: 
The Coulomb matrix elements have been expressed in 



terms of the Racah parameters, and the CEF terms for 
the two-electron case have been determined by those for 
J=5/2. The validity of this j-j coupling picture has 
been carefully dicsussed in the present paper. The /- 
electron hopping amplitudes are obtained using Slater's 
two-center integral. Two types of simplified models have 
also been presented: One is a Hamiltonian including Tg 
and the other for levels. These are believed to be 
useful for further investigations of superconductivity and 
magnetism in /-electron systems 

In this paper the Ts model has been analyzed using 
the exact diagonalization method. We have measured 
several kinds of correlation functions and developed a 
prescription to define the pairing state in multi-orbital 
systems. Furthermore, the symmetry of the orbital de- 
pendent Cooper pair has been carefully discussed. In 
the analysis of the T$ model, it has been suggested that 
unconventional superconductivity should appear around 
the phase boundary between metallic OD-PM and insu- 
lating FO-AFM phases. In addition, we have found evi- 
dence for triplet superconductivity in the ferromagnetic 
phase, induced by strong Hund's rule coupling. With the 
present method of defining the pairing state, it has been 
found that the triplet pairing state is not localized, but 
extended spatially to gain kinetic energy. 

Based on the j-j coupling scheme, it is easy to consider 
the case of several / electrons per site. Although we have 
not shown any results in this paper for the case of n=2, 
corresponding to the U A+ ion (see Fig. 6), we have ob- 
tained that the ground-state is an AFM insulator, which 
is driven by AFM superexchange coupling between neigh- 
boring S=l spins formed by the local Hund's rule cou- 
pling. When we introduce a level splitting into this AFM 
insulating state, we have found robustness for this phase, 
compared with the tijansition around £=0.5 for the case 
of n=l. For UMGa 5 l^fel isostructural with CeMIn 5 , 
if the uranium ions are considered to take U 3+ , again it 
is naively deduced to have an AFM ground state, if the 
local spin is formed by Hund's rule coupling, as shown in 
Fig. 6. However, without any calculations, it is impos- 
sible to determine whether it is a metallic or insulating 
state, although in experiments, metallic AFM phase has 
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been suggested.^! This is left as a future problem. 

As introduced in Sec. 1, PuCoG.su amazingly exhibits 
superconductivity with T c =18.5KJ1j In order to under- 
stand this high T c , first of all, we generally consider that 
the spatial extent of the 5/ wavefunction is large com- 
pared to that of 4/ electrons, leading to a larger hy- 
bridization among conduction and 5/ electrons. Namely, 
the basic energy scale (or effective bandwidth of the 
quasiparticles) in the Pu-compound should be larger than 
that of Ce-materials. Here we envision a scenario simi- 
lar to the high-T c cuprates: If superconductivity in those 
"115" materials is considered to originate from an elec- 
tronic mechanism, T c should be increased as the energy 
scale is enhanced. It may be useful to point out a possible 
analogy with high-T c materials, but we emphasize as a 
separate important issue the active role of orbitals, which 
discussed intensively in this paper. It has been deduced 
from the local-moment behavior that the plutonium va- 
lence is Pu 3+ Jl3 indicating that five /-electrons are ac- 
commodated. In our j-j coupling picture with pseudo- 
spin representation, we can put five electrons in the sextet 
of j=5/2, as shown in Fig. 6, which is just the hole version 
of CeMIns. Thus, the appearance of superconductivity 
itself may be understood without any calculations on the 
basis of our j-j coupling picture. As explained above, 
the drastic changes of T c in the "115" materials may be 
understood, based on the difference in the energy scale 
for 4/ and 5/ electrons, but we believe that another key 



issue related to this problem would be the shape of the 
orbitals. A difference between the electron and hole pic- 
tures should exist stemming from the shape of the active 
orbital, as easily deduced from Fig. 6. If such differences 
among orbitals, i.e., level schemes relating to the CEF ef- 
fect, will be clarified in those materials, it may be possible 
to obtain a unified view to understand unconventional su- 
perconductivity and novel magnetism in "115" materials 
including CeMIns, UMGas, and PuCoGas. This point 
will be further investigated based on the present model. 

As repeatedly mentioned in this paper, there remain 
several open questions in the electronic properties of /- 
electron materials. We believe that the microscopic mod- 
els presented and discussed in this paper will be useful to 
advance our understanding of unconventional supercon- 
ductivity and novel magnetism in /-electron systems. 
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